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Abstract. One of the most intriguing open questions of today's astrophysics is the one concerning the location 
and the mechanisms for the production of MeV, GeV, and TeV gamma-rays in AGN jets. M87 is a privileged 
laboratory for a detailed study of the properties of jets, owing to its proximity, its massive black hole, and its 
conspicuous emission at radio wavelengths and above. We started on November 2009 a monitoring program with 
the e-EVN at 5 GHz, during which two episodes of activity at energy E>100 GeV have occured. We present here 
results of these multi-epoch observations. The inner jet and HST-1 are both detected and resolved in our datasets. 
One of these observations was obtained at the same day of the first high energy flare. A clear change in the proper 
motion velocity of HST-1 is present at the epoch ~ 2005.5. In the time range 2003 - 2005.5 the apparent velocity 
is subluminal, and superluminal (~ 2.7c) after 2005.5. 



1. Introduction 

The giant radio galaxy Messier 87 (M87), also known as 
3C 274 or Virgo A, is one of the best studied radio sources 
and a known 7-ray-emitting AGN. It is located at the 
center of the Virgo cluster of galaxies at a distance = 
16.7 Mpc, corresponding to an angular conversion 1 mas 
= 0.081 pc. The massive black hole at the M87 center has 
an estimated mass — 6 x 10 9 solar masses, with a scale 
of 1 mas = 140 Rs- The bright jet is well resolved in the 
X-ray, optical, and radio wave bands. 

The jet is characterized by many substructures and 
knots. In 1999 HST observations revealed a bright knot 
at about 1" from the core, named HST-1. This feature 
is active in t he radio, optica l , and X-ray regimes. It was 
discussed by IPerlman et al.l (Il999l). who co mpared opti- 
cal and radio images. iBiretta et ail (|l999h measured in 
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Fig. 1. e-EVN image at 5 GHz of the inner jet of M87, 
February 10th, 2010 epoch. The HPBW is 8 x 3.4 mas. 
Levs are -2 2 4 6 8 10 30 50 100 300 500 1000 1500 
mJy/beam 



the range 1994-1998 a subluminal speed = 0.84c for the 
brightest structure (HST-1 East), which appears to emit 
superluminal features moving at 6c. However this motion 
was measured in regions on a larger scale with respect to 
the VLBI structures discussed here. 

VLBI observations of the M87 inner region show a well 
resolved, edge-brightened jet structure. At very high res- 
olution (43 and 86 GHz) near to the brightest region the 
jet has a wide opening angle, and we refer to the many 
published papers which discuss the possible presence of 
a counter-jet and the location of the rad i o core (see e.g . 



Junor et al., 1999; Krichbaum et aUl2005HLv et al .. 2007 



Kovalev et all |2007() 



Very High Energy (VHE) 7-ray emission was reported 
by the High Energy Gam ma-Ray Astronomy (HE GRA) 
collaboration in 1998/99 (|Aharonian et all l2003h . con- 
firmed by the High Energy Stereosco pic System (HESS) 
in 2003- 2006 (lAharonian et~afl 120061) . and by VERITAS 
in 2007 ( Acciari et all 2008 ). Coordinated intensive cam- 
paigns have permitte d to detect the source again in 2008 
(lAcciari et aU l2009h and as recentl y as February and 
April 2010 ()Mariottil I20T0I: lOnd. l2010h . Steady emission at 
MeV/ GeV energies has also been detected by Fermi /L AT 



(lAbdo et all [2009) 



Various models have been proposed to explain the 
multi-wavelength emission and in particular to constrain 
the site of the VHE emission in M87. The inner jet re- 
gion was favo ured by the observed sh ort TeV variabil- 



ity timescales ( Aharonian et al. . 20061) . The VHE emis 



sion could then be produced in the BH magnetosphere 
( Neronov fc Aharonianl | 2007[) or in the slower jet layer 
(jTavecchio fc Ghisellinil . 120081 ). with the spine accounting 
for the emission from the radio to the GeV band; this 
would lead to a complex correlation between the TeV com- 
ponent and the lower energy ones. 
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However, VLBA observations at 1.7 GHz 
20071 ) resolved HST-1 in substruc- 



(jCheung et al 



Aharonian et al 



tures with superluminal components 
(2006) discussed that HST-1 cannot be excluded as a 
source of TeV 7 rays, however they conclude that the 
more promising possibility is that the site of TeV 7-ray 
production is the nucleus of M87 it self. Comparing 
multifrequency data lHarris et al. ( 2008 ) suggested that 
the TeV emission from M87 was originated in HST-1. 

Finally. lAcciari et al.1 ( 20091 ) reported rapid TeV flares 
from M87 in February 2008, associated by an in increase 
of the radio flux from the nucleus, while HST-1 was in a 
low state, thus concluding that the TeV flares originate in 
the core region. 

In this context we started at the end of 2009 a program 
to observe with the e-EVN M87 at 5 GHz to study the 
properties of the M87 core, jet, and HST-1 structure. 



2. Observations and Data Reduction 

Our original monitoring schedule included four epochs 
at 5 GHz, to be carried out before and during the sea- 
son of visibility from the TeV telescopes, namely on 2009 
November 19, 2010 January 2 7, February 10, and March 



2010 ) and April 



28. Followin g the February ([Mariotti . _ 
( One . 2010f) high energy events, three more epochs have 
been added on 2010 March 6, May 25, and June 9 as 
Target of Opportunity (ToO) observations. 

The observations have been carried out in eVLBI 
mode, with data acquired by EVN radio telescopes, di- 
rectly streamed to the central data processor at JIVE, and 
correlated in real-time. The observing frequency of 5 GHz 
was chosen to simultaneously grant a large field of view 
and a high angular resolution. For observations taking ad- 
vantage of the long baselines provided by the Arecibo and 
Shanghai telescopes, our clean beam with uniform weights 
is about 2.0 x 0.9 mas in PA -25°. 

As a result of the large bandwidth (a rate of 1 Gbps 
was sustained by most stations), long exposure (up to 6 
hours per epoch), and extended collecting area, the rms 
noise in our images is mostly dynamic range limited. As 
an average value, we can quote 0.5 - 0.8 mJy beam -1 
in the nuclear region and 0.1 - 0.2 mJy beam -1 in the 
HST-1 region. We present here preliminary results from 
the first 5 epochs. The reduction of May 25th data is still 
in progress. 



3. Results 

3.1. The inner jet region 

In Fig. [T] we show the image obtained on Feb. 10th, 2010 
of the core and inner jet region of M87. The jet orientation 
and velocity has been discussed in many papers compar- 
ing observational data on the j et br ightness and proper 
motion. Recently lAcciari et al. I <l2009h assumed as a likely 
range 9 = 15 - 25 deg. 
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Fig. 2. e-EVN image at 5 GHz of HST-1, January 27th, 
2010 epoch. The HPBW is 8 x 3.4 mas. Levs are -0.3 0.3 
0.5 0.7 1 1.5 2 3 mJy/beam 
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Fig. 3. VLA image at 15 GHz of the inner jet of M87 ob- 
tained on November 2008. The HPBW is 0.10 arcsec. Levs 
are 5 10 15 20 30 50 70 100 300 500 1000 2000 mJy/beam 



We used our images to search for evidence of a pos- 
sible proper motion, comparing different epoch position 
of jet substructures and subtracting images at different 
epochs (with the same grid, angular resolution and simi- 
lar uv-coverage) to look for possible systematic trends. No 
evidence was found in anycase. 

From the jet /counter-jet brightness ratio we derive 
that ficosO > 0.82 which implies 9 < 35°. 

To better constrain the jet velocity and orientation we 
assume that the jet limb-brightened structure is due to a 
velocity structure: a fast inner spine and a relatively slow 
external shear layer. In a range of 9 it is possible that 
the Doppler factor of the slower external regions is higher 
than the Doppler factor of the fast spine. Therefore the 
apparent brightness of the spine can appear fainter than 
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Fig. 4. Distance of HST-1 brightest peak from the M87 core at different epochs 



Table 1. c-EVN Results 



Epoch 


S(core) 


HST-1 peak 


Distance 


yr 


mJy 


mjy/beam 


mas 


2009.87 


1806 


3.5 


905.6 


2010.07 


1810 


2.7 


907.3 


2010.11 


1798 


3.0 


909.5 


2010.17 


1891 


4.6 


911.6 


2010.25 


2013 


3.4 


912.4 



Score uncertainty is ~ 30 mJy; 

Distance is between the HST-1 peak flux and the core, 
with an uncertainty of ~ 2 mas 



that of the shear layer. We estimated different Doppler 
factors using reasonable velocities and jet orientation. We 



found that if the orientation angle 9 is in the range 20° - 
30° the limb-brightened structure can be due to a different 
Doppler factor. 

In Table Q] we report the nuclear flux density measured 
in images at the same angular resolution (8 x 3.4 mas). 
We find a marginal evidence of a nuclear flux density in- 
creasing in the last two epochs, and of a high HST-1 peak 
flux on 2010.17. 



3.2. HST-1 

In our observations HST-1 is clearly resolved (Fig. [2]). 
Comparing different epochs we are able to individuate the 
brightest knot and to measure its proper motion (TableJT]). 
Other features appear to move at about the same velocity, 
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but because of the low brightness and complex morphol- 
ogy we discuss here only the position of the brightest knot 
inside HST-1. 

To better study the dynamic of this structure we 
searched archive VLA data at high resolution (A configu- 
ra tion) and high frequ ency (X, U, and Q bands) . We refer 



to lHarris et a l, (2009) for a discussion of the flux density 
variability. Here we only want to compare different epochs 
to derive the HST-1 dynamic. 

We started t o analy ze data from 2003.6 since, as shown 
( 20091 ). in previous epochs few data are 



Harris et al 



availble, moreover the HST-1 flux density was very low. A 
better analysis with a larger time range will be presented 
in a future paper. 

Starting from 2003.6 the HST-1 structure is well ev- 
ident (see e.g. Fig. [3] obtained on November 2008) and 
well separated by the jet structure near the core. From 
VLA multifrequency observations at the same time we find 
that the total spectrum for HST-1 is moderately steep: on 
2006.21 is afg 4 = 0.68 and a\\ = 0.82 The high frequency 
steepening is not always present, on 2007.52 cv®j 4 is the 



same but 



„15 _ 



0.36. This variability and trend of the ra- 



dio spectrum is in agreement with the substructures and 
variability visible in VLBI images of HST-1. 

We estimated from e-EVN and VLA data the distance 
of HST-1 from the core. In e-EVN data we measured the 
distance between the core and the brightest knot in HST- 
1, in VLA images we used the HST-1 peak, being this 
structure unresol ved. Adding t h e valu es ob tained at 1.5 
and 1 5 GHz by ICheung et all (|2007h and IChang et aT 



( 2010h . respectively, we can study the HST-1 proper mo- 
tion with a good statistic from 2003 to present epoch. 
The apparent proper motion of HST-1 is shown in Fig. 01 
Unfortunatelly in the time range 2006.5 - 2007.5 there are 
only a few points not completely in agreement, however 
taking into account the different observing frequency, an- 
gular resolution and radio telescopes, the proper motion 
looks well defined. 

A clear change in the proper motion velocity is present 
at the epoch ~ 2005.5, coincident with the TeV 7-ray ac- 
tivity and the maximum radio/X-ray flux density of HST- 
1. In the time range 2003 - 2005.5 the apparent velocity 
is 0.5c - 0.6c; in the time range 2005.5 - 2010.25 the ap- 
parent velocity is ~ 2.7c. Assuming a jet orientation angle 
= 25° a proper motion of 2.7c corresponds to an intrinsic 
velocity = 0.94c. 

4. Summary 

With our new e-EVN data obtained in the time range 
2009.87 - 2010.25 we have obtained images of the nuclear 
region of M87 and of the jet substructure HST-1. 

The radio core flux density is constant in the first 
three epochs with an average flux density ~ 1805 mJy 
and slightly increasing in the last two epochs: 2013 mJy 
in 2010.25. 

The inner jet is transversally resolved and assuming 
that the limb-brightened structure is due to a different 



Dopplcr factor because of the jet velocity structure we 
derive a source orientation angle 9 = 20 - 30 degree. 

The HST-1 structure is well resolved in many substruc- 
tures. A proper motion is clearly present. Comparing e- 
EVN data with archive VLA data and published VLBA 
data at 1.7 and 15 GHz we find a strong evidence that 
in 2005.5 HST-1 increased its velocity from an apparent 
velocity = 0.5 - 0.6 c to 2.7c. With present data it is not 
possible to discuss if this change in velocity is related to 
the M87 VHE activity and/or to the maximum radio/X- 
ray flux density of HST-1 at this epoch. A more regular 
and longer monitor and a multi-frequency comparison is 
necessary to clarify this point. 
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